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ABSTRACT 

Finite element, elastic and plastic analyses of various 
aluminum panels, containing holes and notches, were conducted 
for comparison with photoelastic experimental results. A 
FORTRAN IV program, ADINA (Automatic Dynamic Incremental 
Nonlinear Analysis), was used for both linear and nonlinear 
analyses. Mesh refinements were used for each panel and the 
monotonically convergent results were extrapolated using 
Richardson's method. Stresses were locally smoothed from 
the Gauss integration points to the nodal points. Eight 
noded, isoparametric elements were used throughout. Modi- 
Fication to an ADINA preprocessor program, also coded in 
FORTRAN IV, was made for use with a VERSATEC plotter. 

Comparisons were made to the elastic, analytic series 
Bemucion by Howland for a circular hole in a finite strip. 
The finite element results varied by less than one percent 
from Howland's solution. Handbook values for the elastic 
Stress concentration factors of the geometries investigated 
differ from finite element results by less than one percent 
in all cases. The photoelastic works of Frocht were also 
compared where applicable. Stresses in the plastic range 
obtained from slip-line theory for a rigid-perfectly-plastic 
material show excellent correlation to a finite element 
analysis of such a material. Comparisons to elastic and 
plastic experimental data were made for the panels 


analyzed and show good correlation to finite element results. 
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uno UCrT LON 


Development of on-board fatigue monitoring systems 
for Naval aircraft have made it possible to record exten- 
Sive structural loading data in flight. The strain gages 
used in such a system must be located away from stress 
concentration areas to prevent their fatigue; however, 
these areas are of the greatest interest in analyzing 
and predicting fatigue life of the structure. Understand- 
ing the relationship between nominal, far-field stresses 
and local stresses in critical areas thus becomes vitally 
important. Recent experimental investigations into the 
effect of uniform, far-field loads on stress concentration 
areas have been made at the Naval Postgraduate School 
(NPS) using photoelastic techniques [Refs. 1, 2 and 3]. 
These experiments involved loading 7075-T6 aluminum into 
both the elastic and plastic regions, as well as measure- 
ments of residual stresses resulting from plastic yielding. 

Finite element analyses (FEA) of the aluminum panels 
used in the experiments of Stenstrom [Ref. 1] were con- 
ducted. The panels used in the experiments of Engle 
[Ref. 2] and Stuart [Ref. 3] have similar geometry. The 
finite element programs available at NPS were surveyed 
and ADINA [Ref. 4] was chosen for its proven ability to 


produce the nonlinear analyses required for plastic 


5 





yielding of aluminum. To provide increased accuracy, 
each panel was modeled using two meshes. The results 
obtained for the coarse and fine meshes were extrapolated 
to a final result using the Richardson extrapolation 
technique [Refs. 5 and 6]. 

Along with ADINA, a preprocessor program, PSAP1 [Ref. 7], 
was used to verify mesh connectivity prior to analysis by 
ADINA. PSAP1 provides a graphical output of the finite 
element mesh and was coded for the CALCOMP plotter installed 
at NPS prior to 1978. For this thesis PSAP] was adapted 
for use with the VERSATEC plotter now installed at NPS. 

The stress-strain material properties of the 7075-T6 
aluminum actually used to make the panels had to be estab- 
lished to provide an accurate material model for use with 
ADINA. Material testing was conducted to establish the 
Young's Modulus (E), Poisson's Ratio (v), yield stress oe 
strain hardening modulus (E,) and the Ramberg-Osgood co- 
efficients 8 and n. 

Comparisons were made to other works, in addition to 
the experiments conducted at NPS. The initial analysis 
involved a comparison of FEA to the results of Howland 
[Ref. 8], for a circular hole in a finite strip, to vali- 
date the methods used. A comparison of FEA to plane 
stress, Slip-line theory, for rigid-perfectly-plastic 
material was also included as a validation for the plastic 


analyses. 
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Pewee R AGP ROPER Y TESLING@O@F 7075-T6 ALUMINUM 


The elastic and plastic material properties of the 
aluminum panels were established by tensile tests of uni- 
axial specimens made from the same mill run. The speci- 
mens were manufactured and tested according to current 
ASTM standards [Ref. 9]. MICRO-MEASUREMENTS, EA-13-125AD- 
120, precision strain gages with a temperature compensated 
bridge circuit were used on all specimens. Transverse 
gage sensitivity errors were corrected according to the 
manufacturer's recommendations [Ref. 10]. Critical cross- 


section measurements were made with a micrometer. 


Pee LESTS FOR AXIAL LOADING 

Initial tests of the two-gaged specimen, Fig. l, in 
the MTS testing machine indicated a significant bending 
moment was being produced by the 30,000 1b GRIFF grips. 
To investigate this vroblem further, tests were conducted 
on both the MTS and RIEHLE test machines with a five- 
gaged specimen shown in Fig. 2. The results of these 
axial loading tests, shown in Table I, verified that the 
GRIFF grips on the MTS test machine do not give axial 
loading. An inspectionofthe gripped region on the speci- 
men indicated that the jaws of the grip were not applying 


a uniformly distributed force and thereby induced a bending 


iS 





moment by off-axis loading as shown in Fig. 2. The grips 
on the RIEHLE test machine gripped evenly and a uniform 


strain distribution resulted as seen in Table I. 


B. CHARACTERISTICS OF 7075-T6 ALUMINUM PANELS 

The following characteristic properties were determined 
from the four specimens tested. 

fm round Ss Modulus (5) 

Tests were conducted using the specimen shown in 
fies) Ss On the MTS test machine with 10,000 lb INSTRON grips, 
which gripped the specimen evenly. The results of testing 
three specimens are shown in Tables II to IV. Linear 
regression in the elastic range of all the test results 
determined a Young's Modulus of 10.12 x 10° psi, with a 
Serrelation coefficient of 0.9996. 

2. Poisson's Ratio (\v) 

Tests were conducted using the specimen shown in 
meg. | On the RIEHLE test machine with 10,000 1b RIEHLE 
grips. The results are tabulated in Table V. Linear 
regression of these results in the elastic region deter- 
Mined Poisson's Ratio to be 0.3256 with a correlation 
Seerficient of 0.99996. 

3. Yield Stress and Strain Hardening Modulus 

These values, required for ADINA's bi-linear 

material model, were determined graphically using the 


data from Tables III and IV. Plastic region data in 
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Table II is not reliable because of excessive creep en- 


countered during that test. 


0.2% offset yield stress, o,, 76,000 psi 


roo, UG On sd 


strain hardening modulus, EY 
The graphical fit of these values to the test data can be 
seen in Fig. 4. 

4. Ramberg-Osgood Coefficients 


The Ramberg-Osgood equation for elastic-plastic 


Stress-strain characterization 1s given by: 


Oo o,n 
= — + — 
e = 8 (5) (1) 
where: 
ge = strain 
Oo = stress 


E = Young's modulus. 
The 8 and n coefficients were determined graphically from 
the data of Table IV, by the method given by Rivello [Ref. 11]. 
The data in Table IV gave the following values which are 
the best fit to the combined test data 


eo dae 


n 24. 58 
Mees graphical fit of these values, in Eq. (1), with the 


test data is also shown in Fig. 4. 
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ieee emoue LEAT TON TO GRAPHICAL PREPROCESSOR 


The use of a graphical preprocessor program, such as 
PSAP1, is vital in detecting mesh errors that may other- 
wise go unnoticed. Establishing the proper node locations 
and element geometry prior to analysis for a complex code 


such as ADINA is of utmost importance. 


A. PSAP1L MODIFICATIONS 

The program PSAP1 was originally coded in FORTRAN IV 
for use on the NPS IBM 360/370 installation with the 
CALCOMP Model 765 drum plotter. The CALCOMP system once 
installed at NPS used the +Y axis as the unlimited plotting 
axis, see Fig. 5. The entire plotting logic in PSAPI1 uses 
this orientation of axes to allow multiple plots in a con- 
tinuous strip. With the VERSATEC Model 8222A electrostatic 
plotter now installed at NPS, the +X axis becomes the un- 
imemieced plotting axis, shown in Fig. 5. To avoid an ex- 
tensive recoding of PSAP1 for use with the VERSATEC plotter, 
a Simple coordinate transformation of the plot was made in 
a limited number of short subroutines. First, all installa- 
tion dependent plotting calls used in PSAP1 were identified. 
These involved seven plotter functions for which new sub- 


melitines were coded. 
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rumet lon New Subroutine 


Initialize Plotter CALCMP 
Move Plotter Pen CALPLT 
Letter on Plot NOTATE 
Number on Plot CALNUM 
Determine Current Pen Location CALWH 

Draw a Line CALINE 
Stop Plotter PSTOP 


The subroutines listed above merely rotates the plot 
to coincide with the VERSATEC axis orientation and retain 
all features originally in PSAP1l. Since all plotter 
hardware code is now isolated in these seven subroutines, 
future adaptations to other plotting systems is simplified. 
To provide documentation of this update to PSAP1, a complete 


listing of the new program is provided in Appendix D. 


pee UoE OF PSAPIL 

Previous use of PSAP] on the IBM 360 system necessitated 
use of a load module since PSAP1 took over one minute of 
PU Dine to compile. With the new IBM 3033 system 
compilation requires eight seconds; however, use of a load 
module or disk stored source code is still recommended 
Since PSAP1 contains roughly 2,500 lines of code. Appendix 
A contains sample JCL to use PSAP1 on the IBM 3033 MVS 


system. With minimal effort PSAP1 could also be set up 


Ike, 





for use on the IBM VM/370 system. The user's manual for 
PSAP] is in Ref. 7. In addition to a mesh plot, PSAP1l 
provides a listing of the node coordinates, element 
connectivity and several key input values used in execu- 


tion of ADINA. This information provides a useful check 


@f the input data. 
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Lael Te ELE MENGANALYSIS (PEA) 


free DESCRIPTION OF MESHES USED 

1. Length to Width Ratios and Boundary Conditions 

Initial finite element models of specimens had 

length to width ratios near one, as used by Garske [Ref. 12] 
for his FEA, but they did not provide the desired uniform 
distribution at the loading boundary. Specimen length to 
width ratios of 3-5 were used by Armen, Pifko and Levine 
[Ref. 13] in their FEA and by Stenstrom [Ref. 1] in his 
photoelastic experiments. The criteria established to 
determine uniform boundary stress distribution was uni- 
formity in nodal displacements along the loaded edge as 
discussed by Segerlind [Ref. 14]. In the models used for 
FEA in this thesis, nodal displacements were uniform to 
Within 0.1%, and the resulting stress distribution was 
uniform axially to within 0.1% at the panel ends. [In all 
cases two-dimensional, eight noded, isoparametric elements 
were used. These higher order elements cannot be loaded 
in an "intuitive" manner as discussed by Zienkiewicz 
Meet. 15, p. 223]. Figure 6 shows the nodal loading 


required to obtain a uniform surface load. 
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2. Element Meshes 

Two meshes were developed for each panel analyzed. 
A reasonable effort was made to keep element corner angles 
as close to 90° as possible to reduce the effect of element 
distortion discussed by Hopkins and Gifford [Ref. 16]. 
All meshes modeled a guarter of the actual panel by using 
the two axes of symmetry aS 1S common practice in FEA. 
The step from course to fine element meshes was made so 
that each element in the course mesh was subdivided into 
four smaller elements of the same type. Such a mesh sub- 
division can be expected to give monotonic convergence of 
results, Cook [Ref. 17], and allow extrapolation to results 
Sean infinitely fine mesh. Figures 7 through 13 illus- 
trate the element meshes used in this analysis as plotted 


by PSAP1. 


Eee COMPUTATIONAL PROCEDURES 
1. Using ADINA 

Once the mesh has been developed, input data is 
prepared in accordance with the ADINA user's manual [Ref. 4]. 
This same set of data is then used as input for PSAPI1 to 
check for errors and provide a graphical display of the 
element mesh. After preprocessing by PSAP1, the data is 
entered into ADINA for analysis. In the case of linear 
analysis, two types of stress output may be specified, 


nodal point or Gauss integration point. Nodal point output 


a 





can be computed for up to eight node point stresses for 
each element. Since 2x2 Gauss integration was used, four 
Gauss point stresses were computed for each element. The 
2x2 Gauss integration is recognized as the most efficient 
integration order for this type of analysis [{Ref. 15, p. 284]. 
The linear analysis used an isotropic linear elastic material 
model (MODEL "1" in Ref. 4) which required input of E and 
v Material properties. The nonlinear analysis allows only 
Gauss point stress outputs and uses a bilinear elastic- 
plastic material model, with von Mises yield condition and 
isotropic strain hardening (MODEL "8" in Ref. 4). 

For static analyses ADINA uses a time function 
method to apply loads in steps. Linear analysis loading 
was accomplished in a single step to a nominal value of 
3,000 lbs load. Nonlinear analysis loads were applied in 
ten steps to a maximum value, matching the experimental 
loads, and then unloaded to zero in ten steps to obtain 
residual stresses. The stress output from ADINA is a 
listing of nodal or Gauss point stresses for each element. 
Since the only area of interest in this analysis was the 
distribution of stresses along the reduced cross-section, 
no large post-processing program was developed or used. 
All final computations using ADINA output data were accom- 
plished on a HEWLETT-PACKARD 9830A calculator, using short 


programs coded in BASIC. If more extensive stress 
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distribution information were desired, some form of auto- 
mated Re ereeicessing would be necessary to reduce the 
computational workload. At a minimum, nodal stress outputs 
by ADINA must be averaged to obtain unigue values of stress 
at nodes shared by more than one element. 
2. Richardson Extrapolation 

The use of course and fine meshes allows extrapola- 

tion to an infinitely fine mesh as discussed earlier. 


Richardson extrapolation [Ref. 5] was used in this analysis 


where: 
m ™m 
: o, (AS) oh) ing 
= eGap m m 
h -h 
r Se 
where 
Fextrap = extrapolated solution 
oe = solution obtained with h=h | 
Oe = solution obtained with h=h, 
A. = linear dimension of course element 
he = linear dimension of fine element 
m = 2 (for this analysis) 
The exponent m is determined by the order of the discre- 
tization error O(h'). Since h represents the length of an 
element the element area is represented ae In a two 


dimensional problem such as this O(h") is of the order of 


h*, the area of an element. In the mesh refinement scheme 
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h 
ae Ey Selle . ; 
used ne = 7h or h Sor PeuUatlen s(2 scan be rewritten 
h i 
EF, 2 c, 2 
Oo i = oP A? 
Yextrap ~ eer ee ey 
Sa 
h h 
S E 
thus hes l : 
aT 
5 aes ee. c (4) 
Speco 3 
- 


Equation (4) then becomes the relation to obtain extrapolated 
stresses from coarse and fine mesh results in a two dimen- 
Sional analysis. Better extrapolations can be obtained by 
uSing three or more refined meshes, but the computational 
effort increases Significantly. 

3. Optimal Stress Locations and Local Smoothing 

It 1S generally accepted that the most accurate 

Sampling points for stresses are the Gauss integration 
points within the element [Ref. 15, p. 281, and Ref 18]. 
In this analysis, the nodal points are of the greatest 
interest; thus a technique of local smoothing must be 
applied to the integration point stresses to obtain nodal 
stresses as reported by Hinton and Campbell [Ref. 19]. 
The formulation of this local smoothing technique for ADINA 
elements is developed in Appendix B. The nodal values 
obtained must then be averaged if shared by two or more 


elements. 
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4. Computational Times 

Because of the extremely large size of ADINA (about 
17,000 lines of code in the NPS version) and the out of core 
solver, it does not adapt well to time sharing systems. 
Using the IBM 360 system at NPS, ADINA required 31 user 
defined overlays to create a manageable load module in 
about 30 minutes of CPU time. With the new IBM 3033 MVS 
system at NPS, ADINA is compiled without overlays in about 
one minute. When using a load module, the program execu- 
jmren took less than 2 minutes CPU time on the IBM 3033. 

In addition to ADINA, the preprocessor (PSAP1) 
and post-procesSing techniques involve considerable time 
eaemeertrort. Figure 14 is a flow chart of the computational 
procedure used in this analysis. An example of the JCL 
to use ADINA in load module form on the IBM 3033 MVS 
svstem with use of the mass storage facilities is shown 


maeADppendix C. 
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VY. RESULTS OF ANALYSIS 


fee)6C LRCULAR HOLES IN LINEAR MATERIAL 

The FEA results for a circular hole in a finite width 
Strip were used to validate the elastic computational pro- 
cedure discussed earlier. The results of Howland [Ref. 3] 
were compared to both the Gauss point smoothed results and 
the nodal output resultS in Fig. 15. The stress concentra- 
@ug@emractOr 0/O— is rererenced to the far-field stress. 
The smoothed results give the best match to the results 
of Howland at the edge of the hole, and the only sSignifi- 
Cant variation between the two FEA methods occurs within 
the first 0.25 inches from the edge. In order to obtain 
the 0.25 inch stress value for the coarse mesh, in the 
Gauss point smoothed result, a midside node value had to 
be obtained by the averaging method discussed in Appendix 
B. The linear distribution of smoothed stresses alcng the 
Sides of the element, [Ref. 19], appears to produce a 
less accurate result in this area of extreme stress 
gradient, when compared to ADINA's nodal output result. 
This tendency was noted in all cases; however, the peak 
stress values from smoothed results consistently gave 
better correlation with other investigators [Ref. 20]. 

A circular hole with \A=0.25 was also analyzed and 


compared to the experimental results of Stenstrom [Ref. 1] 


2a 





along with an interpolation of Howland's results. The Oy 
experimental data correlates well with the FEA results; 
however, the os experimental data shows significant varia- 


tion between 0.125 and 0.375 inches from the edge of the 


hole, as seen in Fig. 16. 


Bee OPPOSITE U NOTCHES IN LINEAR MATERIAL 
The results for linear analysis are presented in non- 
dimensional stress concentration form; however, the 
normalizing stress changes. For U notches Km 1s the 
stress concentration factor referenced to the theoretical, 
nominal stress (Oo) in the reduced cross-section where 
= Load/Area of Reduced Cross-Section. 
ioe ohnallow Notch Panel 
The FEA results plotted with the experimental 
Mata Of Stenstrom are shown in Fig. 17. Once again both 
FEA results are shown and the variation for the two methods 
occur within 0.25 inches from the notch edge; however, 
there was less variation than was seen in the circular 
hole analyses. 
For this panel the experimental data appear to be 


uniformly below the FEA results for o The O,. data shows 


5° 
Seemicicant variation at the 0.625 inch point but follows 
the proper trend within 0.5 inches from the notch edge. 


According to data collected by Peterson [Ref. 20], this 


notch geometry should yield a maximum Kan = 2.74. The 
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GaussS point smoothed results matched this value exactly. 


Stuart [Ref. 3] reported a as 2.69 for this same notch 


geometry, with a standard deviation of 0.187 in the 14 
samples he measured by photoelastic methods. Early pnoto- 
elastic work by Frocht [Ref. 21] determined Ke = 27 ~£or 
this notch geometry. The FEA results appear to be in good 
agreement with other investigators, for this notch geometry. 
2. Deep Notch Panel 

Results of the FEA for a deep notch were plotted 
with Stenstrom's experimental data in Fig. 18. The results 
of the two FEA methods again diverge within 0.5 inches from 
the notch edge. FEA stress values at the 0.25 inch point 
have spread farther apart in this case Since the stress 


gradient iS very severe at that point. The experimental 


data correlates well for both Ta and Oi however, the maxi- 


mum experimental o, is considerably lower with a Kan = 95 903 


0 


Results reported by Stuart for this notch geometry was a 


Kr = 94.05 with a standard deviationof 0.219 for 14 speci- 


mens measured by photoelastic methods. Frocht [Ref. 22] 


reported a photoelastic Kn wom TrOlesHlS notch geometry, 


but concluded that the result was 5-10% low, giving a 


corrected range of Kan from 4.1 to 4.3. The Gauss smoothed 


FEA result gave a K, = 4.24 which compares well to an 


a 
empirical relation given by Peterson [Ref. 20] for 


mar< 0.25. 


ae 





K=(1 - 5) (0.78 + 2.243vE72) 0.99: + 0.18 <= - 1.06 (24) ° + 1.712%? | (5) 
where 
me= riotch aepth (3.9375) 
fe Oren radius (0.625) 
ad = minimum width (15.625) 
D = maximum width (23.5) 


Inserting the values above into Eq. (5) gives a K, = 4.26, 


T 
which is 1/2% above the FEA result. Other FEA results 
reported by Armen, Pifko and Levin [Ref. 13] and Griffis 
[Ref. 23], using linear strain triangle (LST) elements, 


produced K,, values within 5% of those produced by use of 


oe 
Eq. (5). The rectangular elements used in this analysis 
are known to give better results than LST elements as 
noted by Clough [Ref. 24]. It is clear that the FEA 


results obtained for this notch are in good agreement 


With other works. 


C. OPPOSITE U NOTCHES IN NONLINEAR MATERIAL 

The analysis for loading into the plastic region of the 
7075-T6 aluminum was made uSing the bilinear material 
model discussed earlier. The loads used were selected to 
match those used in the experiments of Stenstrom; thus 
allowing direct comparison. The strains obtained in those 
experiments were used to solve for stresses by use of the 


Prandtl-Reuss plastic flow equations. 
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1. Shallow Notch Panel 
The results of FEA for the three load cases, 60,000, 
65,000 and 70,000 lbs are presented along with the experi- 
mental results in Figs. 19 through 21. The Op results 


compare well although no trend for peak o, stress away from 


9 
the notch edge is shown in the experimental data. In all 


cases the FEA determined the peak o, stress to occur near 


8 
the yield boundary, and the gradient of the o, stress to 


8 
mee OLf dramatically in the plastic zone. This character- 
istic behavior of the d, stress was reported by other 
investigators [Refs. 13 and 23] using FEA on 2024-T3 aluminum. 
mame elastic-plastic stress distributions reported by 
Frocht [Ref. 25] show similar trends. The experimental 
data also shows a marked change in the gradient ofc, stress 
within the plastic region. The growth of this plastic 
region iS approximated using the FEA results for this notch 
meeeigs. 22 through 24. 

Experimental data for the Oo, stress Gasset don 

matches the FEA results closely except at the notch edge 
where the measured O.. does not go to zero as it should. 
The characteristic peak value of Oo. near the plastic 
boundary as seen in the FEA results is also shown by the 
Eest data. 

The FEA residual stress computed upon unloading 


From the three load cases are shown in Figs. 25 through 30. 


Mio characteristic distributions of the 4 residual stress 
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agrees with those reported by others [Refs. 25, 26 and 27]. 
The experimental residual stress distributions reported by 
Stenstrom show similar trends but significant variations 
when compared to the FEA results. 
2. Deep Notch Panel 

Three load cases were computed to plastic loading 
levels; however, only limited experimental results are 
available for this notch as seen in Fig. 3l. The 30,000 1b 
load is just at the onset of yield in the notch root area. 
Limited residual experimental data [Ref. 3] was available 
for comparison in Figs. 32 and 33 which are plots of the 
residual Op and o stress distributions as a result of the 
three loading cases. The 30,000 1b load has caused yielding 
in a small region at the root of the notch as seen in Fig. 34. 
Figures 35 and 36 illustrate that the plastic zone does not 
grow to the extent it did in the shallow notch. The stress 
gradients are very severe in the deep notch and as a result 
the 0.25 inch sampling points used in the FEA may be useful 
in only showing gross trends close to the notch edge. The 
trends appear to be much the same as in the shallow notch, 


with the peak g. and oe stresses occurring near the yield 


9 
boundary. The yield boundaries shown in Fig. 31 are 
approximations based on qualitative analysis of the finite 


element results. The - experimental data for the 30,000 


lb load case correlates well with the peak stress, again 
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appearing low as it did in the linear analysis. The residual 


o, stress distribution in Fig. 32 follow much the same trends 


9 
as seen in the shallow notch but with much higher gradients 
within the first 0.5 inches from the notch. The FEA results 
for the residual co, stresses shown in Fig. 33 indicate a 
limitation in the element size used since the Oe value of 
the notch edge does not return to zero as it should. Because 
of this problem the data may be questionable for showing 
proper trends in the first 0.5 inch from the notch edge. 
3. Rigid-Perfectly Plastic Panel 

A stress distribution for the theoretical material 
used in slip-line theory was desired. By using ADINA's 
bilinear material model such a material could be approxi- 
mated reasonably well. For this analysis a Young's Modulus 


216 


my Of 10 pSi was used to model perfect rigidity. fThe 


Strain hardening modulus (E,) was set to zero to model 


ie 
Merrect plasticity. Poisson's ratio (v) was 0.4999999, 

as close to 0.5 as the computer would allow. The E = ome 
also represents a computer limit in approaching an infinitely 
large E. Figure 37 illustrates the results obtained and 
compares them to a slip-line solution. The results are 
normalized to the arbitrary 73,000 psi yield stress used 

in the analysis. The T , values obtained agree exactly 


with slip-line theory. Conversely the C,. results do not 


reflect the same values as slip-line theory, but do show 


S35 





a similar trend. The growth of the plastic zone obtained 


ms shown in Figs. 38 through 40. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The results obtained from the FEA have proven useful 
in determining the validity of experimental data gathered 
by photoelastic techniques. The FEA results varied by 
less than 1% when compared to published analytical results 
and handbook values. Experimental data from Stenstrom's 
photoelastic work correlated well with the FEA results. 
The primary exception was the residual stress experimental 
values, which varied significantly from the FEA results. 
The variation may be in transforming the residual strains 
measured photoelastically into stresses for comparison 
with the FEA results, since ADINA only provides a stress 
output. Limitations of ADINA's bilinear material model, 
initially considered severe, do not appear to have hampered 
this investigation. A possible exception is the residual 
analyses where the transition region from elastic to 
plastic strains becomes especially important. The Gauss 
point smoothed results gave the best correlations at the 
edge singularities in all cases; however, due to the 
limitations noted at 0.25 inches from the edge, the 
results at that specific point may not be as accurate for 
this method. The nodal output results gave consistently 


higher stress values at the edge singularities. Because 
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of the severe stress gradients near the deep notch analyzed 
the use of a finer element mesh near the notch would probably 
produce better results. 

The effort involved in developing two meshes such as 
those used in this thesis is considerable. An automatic 
mesh generation capability would reduce the workload and 
allow experimentation with several types of element meshes. 

ADINA proved to be a useful and powerful program, as 
expected, but something simpler and less awkward to use may 
be all that is required for two dimensional analysis. Such 
a system is already in use at NPS but does not offer non- 
linear capabilities. If use of ADINA is to be continued in 
this type of investigation, a post-processing program should 
be adapted. There are programs available to post-process 
ADINA data at NPS [Refs. 28 and 29] but they would require 
modifications to work with two dimensional analyses and 
the VERSATEC plotter. 

Standardized material property testing would ease the 
inevitable task of obtaining basic material properties for 
use in analysis or experiments. Some form of automatic 
data collection with use of the MTS testing machine would 
allow testing of a larger sample population and provide 


statistically more accurate information. 
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EXAMPLE OF COMPLETE PANEL MESHES 
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TABLE I. MTS AND REIHLE 5 GAGE TEST RESULTS 
: ee 
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TOBE sll. MTS SPECIMEN A TEST RESULTS 
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TABLE Lit. MrS SPECIMEN B TEST RESULTS 
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TABLE IV. MoS, SeECUMEN, C TEST RESULTS 


Cross~section = 0.03975 in 
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TABLE V. RETHLE SPECIMEN TEST RESULTS 
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TABLE VI. X = 0.2 HOWLAND DATA 


DISTANCE 
FROM HOLE 87 36 





TABLE VEL: X = 0.25 HOWLAND DATA 
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Aes tore eee. A = 0.2 FEA RESULTS - NODAL OUTPUT 
Cte 5000 psi 
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TABLE IX. hee= eee HEA RESULTS = GAUSS OUTPUT 
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RABBLE Xx. hes 0.25 FEA RESULTS = NODAL OUTPUT 
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TABLE XIII. SHALLOW NOTCH FEA LINEAR RESULTS - GAUSS 
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TABLE XIV. DEEP NOTCH FEA LINEAR RESULTS - NODAL 
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TABLE XVI. SHALLOW NOTCH FEA NONLINEAR 60,000 1b LOAD 


DISTANCE NO LOAD RESIDUALS 


FROM NOTCH,in. oy oss S27, DSi Oy ;PS1 - ,PS1 


HH, SHS 4 PF liad | § weet | 1, oh dts, lf 


-_ -_ ‘= -— —_— “-_ 


ye a Meee F523, 5 <b SS ~F6So, | 
a. 50 eee Lease. HI DLSe, 
aH. 75 Shs4e, 2 Lsa03 4 2446. aia 
1.80 Se4i1.4 LSTsa.¢ 1age.8 tPa6.7 
1.25 S2STH. 1 Sid7. fier, iSGi.é 
1.50 49111.1 17477. 1 a25,7 1258.3 
1.75 40088.7 te65e7.2 Sad,4 | ta4ae.7 
a0 42573, 3 iseer.s Fda a25. 0 


Pasir XVII. SHALLOW NOTCH FEA NONLINEAR 65,000 1b LOAD 


NO LOAD RESIDUALS 


DISTANCE 
FROM NOTCH,in. oy ;pSL Oo ,psi te ,pSi or »psi 








TABLE XVIII. SHALLOW NOTCH FEA NONLINEAR 70,0001Lb LOAD 


NO LOAD RESIDUALS 





DISTANCE 


FROM NOTCH ,in. Tp ;PpSl oie ,PS1L Tp POS: oe ps1 





_f 
fut 
i} 

J 
74 
= 
4 ji 
Ee 

oF 
iT; 
Naa 

- 
ae 
it 
= 
ji 
i at 
[ua 
i 
=] 
ma 


TABLE XIX. DEEP NOTCH FEA NONLINEAR 30,000 1b LOAD 
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MABLE XXII. REGED -  PERPEecirby = PLASTIC RESULTS 
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TABLE XXIV. EXPERIMENTAL DATA SHALLOW NOTCH LINEAR LOADING 
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TABLE XXV. EXPERIMENTAL DATA SHALLOW NOTCH 60,000 lb LOAD 
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TABLE XXVI. EXPERIMENTAL DATA SHALLOW NOTCH 65,000 lb LOAD 
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APPENDIX B 


Ss LOCAL LEAST SQUARES SMOOTHING 


® corner Nodes 
o Midside Nodes 


+ 2x2 Gauss Points 





Two-Dimensional Isoparametric Element from ADINA [Ref. 4] 


The local smoothing expression from Hinton and Campbell 


[Ref. 19] in ADINA coordinates becomes 
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The midside node stress values may be obtained by averaging 
the values at the associated corner nodes, since the distri- 
bution of the smoothed stresses is linear along the sides of 
the element. Smoothed stress values obtained by this least 
Squares method should subsequently be averaged to obtain 


unique values at nodal points shared by adjacent elements. 
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